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3-D  DISPLAY  SYSTEM  ANALYSIS  AND  TEST 
ABSTRACT 

This  Teclmical  Report  presents  the  results  of  various  analyses  and  tests 
conducted  at  the  Naval  Training  Device  Center  during  the  summer  of  1970 
concerning  the  3-D  Display  System  supplied  by  the  Bowla*  Fluidics 
Corporation. 

The  nonlinear  differential  equation  of  flow  from  a  fluidic  pulse  jet  motor 
(FPJM)  is  presented  and  solved  for  the  steady  state  periodic  jet  velocity 
caused  by  a  periodic  pressure  variation.  An  expression  for  average  FPJM 
thrust  is  given  as  well  as  a  WANG  700  computer  program  for  actual  thrust 
computation. 

The  depth  stability  of  a  submerged  model  using  the  Cartesian  diver 
principle  was  investigated  and  the  differential  equation  of  motion 
developed.  A  depth  control  system  was  devised  and  is  described. 

Thrust  measurements  were  conducted  on  models  using  the  Bowles  supplied 
hardware.  The  measurements  Are  in  good  agreement  with  the  analytical 
results  produced  by  the  WANG  700  computer  program. 
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SECTION  I 


INI  INDUCTION 

Initial  development  of  a  unique  three-dimensional  display  system, 
figure  1,  was  reported  in  Bowles  Fluidics  Corporation  report.  Reference  1, 
and  hardware  was  delivered  to  the  Naval  Training  Device  Center,  Orlando, 
Florida  in  May,  1970  undet  Contract  N61339-69-C-0078. 


Figure  1.  Tank  and  Model 

The  system  consists  of  a  water  filled  tank  containing  a  submarine  model 
which  moves  in  three  dimensions  in  response  to  water  pressure  variations 
caused  by  motion  of  a  hydraulically  driven  piston  in  a  cylinder  located  in 
the  side  of  the  tank,  figure  1. 
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Pigure  2.  Tank  Cylinder  &  Piston  Hydraulic  Actuator 


Motion  of  the  piston,  in  turn,  is  controlled  by  a  feedback  system  with 
three  inputs:  (1)  Average  position,  (2)  Amplitude  and  (3)  Frequency  of 
periodic  motion,  figure  3, 


2 


Figure  3.  Input  Controls 


(a)  On-Off  Switch 

(b)  Average  Piston  Poation  (Depth  Control) 

(c)  Amplitude  of  Platon  Motion  (Velocity  Control) 

(d)  Frequency  of  Piston  Motion  (Turning  Control) 


Vertical  position  of  the  model  Is  controlled  by  the  average  tank 
pressure  using  the  principle  of  the  Cartesian  diver  (devil  or  laq>).  A 
Cartesian  diver  is  a  toy  consisting  of  a  small  figure  sealed  in  a  liquid 
filled  tube  with  vertical  position  controlled  by  finger  pressure  on  s 
diaphragm  or  cork  sealing  the  top  of  the  tube. 

Horizontal  velocity  of  the  model  is  controlled  by  the  amplitude  of  a 
periodic  tank  pressure  change,  and  model  turning  to  the  right  or  left  is 
controlled  by  the  frequency  of  the  periodic  pressure  charge.  Model  turning 
is  possible  because  the  discharge  tube  axes,  see  figures  A  and  5,  do  not 
pass  through  the  model  center  of  gravity  and  any  unbalance  in  flow  from  the 
two  tubes  causes  an  unbalanced  jet  thrust  which,  In  turn,  provides  a  turning 
jjoment  on  the  model. 
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Figure  5.  Display  Model  (Side  View) 
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Control  of  the  relative  flow  from  the  two  discharge  tubes  is  obtained 
by  adjusting  the  frequency  of  the  tank  pressure  variations  over  a  range 
from  the  tesonant  frequency  of  one  tube  to  the  other.  For  the  model  s'aown 
in  figure  4  the  maximum  right  turning  moment  occurs  when  the  tank  pressure 
is  varied  at  8  cycles/sec  ana  maximum  left  turning  moment  occurs  at 
11  cycles/sec. 

The  different  resonant  frequency  for  nozzle  flow  from  the  two  sides 
depends  upon  a  stiff  vertical  longitudinal  partition  which  divides  the 
model  into  a  port  and  starboard  side.  Each  side  contains  a  fixed  trapped 
air  volume  of  0.4  cubic  inches  which  resonates  with  the  fluid  flowing  in 
the  corresponding  exit  tube.  The  tubes  are  of  sufficiently  different 
lengths  to  provide  the  observed  difference  in  resonant  frequencies.  Figure 
6  shows  model  construction  details,  and  an  analysis  of  the  thrust 
produced  by  these  pulse  Jet  motors  is  given  in  Section  IV. 

The  system  allows  the  model  to  be  maneuvered  throughout  the  tank  with 
no  visible  control  devices  such  as  vires  or  strings  and  with  the  desired 
motion  determined  solely  by  means  of  fluid  pressure  variations. 


JCjDx  A&L  BAASS  "TUBE 

32T 


Figure  fa.  Sketch  of  Submarine  Model 
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SECTION  II 

STATEMENT  OF  THE  PROBLEM 

Control  of  the  display  system  model  by  one  operator  Is  difficult 
because  the  depth  control  Is  unstable.  An  investigation  was,  therefore, 
undertaken  to  establish  design  criteria  for  a  suitable  depth  control 
stabilizing  system.  It  was  quickly  discovered  that  the  thrust  performance 
analysis  for  fluidic  pulse  Jet  motors  as  presented  In  the  Bowles  report, 
reference  1,  needed  to  be  expanded  before  it  could  serve  as  a  basis  for 
a  detailed  stability  analysis.  This  report  presents  the  results  of  the 
thrust  and  depth  stability  analyses.  Experimental  results  are  given 
as  well  as  suggestions  for  future  effort.  A  computer  program  for  the 
Wang  700A  calculator  which  may  be  used  to  compute  the  thrust  of  s 
fluidic  pulse  jet  motor  is  included  as  an  appendix. 
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SECTION  III 
LETTER  SYMBOLS 

A  ■  nozzle  tube  flow  area,  in^. 

B  -  viscous  damping  coefficient,  lbf-aec/ir. 

C  *  capacitance,  farads  (or  in^/lbf.)* 

CD  •  drag  coefficient,  dimensionless 

d  *  differential  operator 

D  or  Dn  ■  nozzle  tube  diameter,  in. 

Dd  ■  effective  diameter  of  diaphragm,  in. 

e  ■  electromotive  to rce,  volta  (or  Ibf/in^.)* 

f  ■  frequency,  cycles/sec.  ,  or  hertz 

F  ■  thrust,  Ibf. 

»  buoyant  force,  Ibf, 

F  ■  force  of  gravity,  Ibf. 

© 

Fffl  *  motor  thrust,  Ibf. 

F^j,  ■  motor  thrust  with  no  nozzle  throttling,  Ibf. 

Ffl  -  shear  force  due  to  viscosity,  Ibf. 
g  ■  acceleration  of  gravity,  in/sec. 

gc  *  constant,  386  (numerically  equal  to  standard  acceleration  of 
gravity  in  units  of  in/ sec*) 

i  •  current,  amperes  (or  in/sec.)* 

j  -  FT 

1  *  nozzle  tube  length,  in. 

2 

L  •  inductance,  henrys  (or  Ibm/in  .)* 
n  ”  exponent 
N  -  R/|U[ 

m  *  mass,  lbc. 

o  ■  resonant  condition  -  uced  as  subscript 
♦Units  for  analogous  mechanical  quantities 
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P  ■  pressure,  lbf/in2. 

2 

Pn  ■  ambient  pressure  external  to  model,  lbf/in  . 
P  *  pressure  phasor,  lbf/in2. 


q 
Q 

R 
S 
t 
T 

u 
U 

V  or 

Vs 

x 

X 


dynamic  pressure , /l/2)  —  u2 ,  lbf/in2. 

2* 

selectivity 

2 

resistance,  ohms  (or  lbm/in  .  -  sec.)* 

2 

surface  area  of  nozzle  interior,  in  . 
time,  sec. 
period,  sec. 

Jet  velocity,  in/sec. 
velocity  phasor,  in/sec. 

3 

Va  ■  trapped  air  volume,  in  . 

3 

volume  of  solid  parts  of  diver,  in  . 
nozzle  opening,  in. 

average  nozzle  opening  over  a  period,  in. 
vertical  dimension,  in. 


<X 

fi 

A 

$ 

s 

u) 


/rnf'Dm 

incremental  operator 
incremental  operator 

3 

fluid  density,  lbm/  in  . 
Ziff‘s  angular  velocity 
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SECTION  IV 


FLUIDIC  PULSE  JET  MOTOR 

A  fluidic  pulse  jet  motor  (FPJM)  Is  represented  by  figure  7,  It 
consists  of  a  chamber  with  a  trapped  air  ?pace  of  volume  V  and  a  jet  nozzle, 
shown  for  simplicity  as  a  symmetrical  reentrant  tube,  of  inner  diameter  D 
and  length  j[  .  When  the  surrounding  ambient  fluid  pressure  is  varied 
periodically,  fluid  alternately  enters  and  leaves  the  nozzle.  The  rate  of 
change  of  momentum  imparted  to  the  fluid  causes  a  reaction  force  to  be 
developed  on  the  motor.  No  steady-state  rate  of  change  in  momentum  is 
imparted  to  entering  fluid  since  whatever  axial  velocity  is  present  in  the 
nozzle  is  destroyed  upon  impact  with  the  chamber  inner  surface.  Thus,  the 
time  average  motor  thrust  may  be  shown  to  be  proportional  to  the  average 
momentum  of  the  fluid  leaving  the  motor  per  unit  time.  In  this  respect 
the  motor  is  the  analog  of  a  rectifier  in  an  ac  electrical  circuit. 


chamber 


AIR  VOLUME,  V 


4  ambient  pressure. 


—  Fluid  Pressure, 
Density.^* 


-  -  -  - - —  Jet  velocity,  /u. • 


Jet  nozzle 
area,  A 
diameter,  D 

length,  X 


Figure  7.  Fluidic  Pulse  Jet  Motor 


4.1  Differential  Equation 

An  analysis  of  a  FPJM  can  be  made  by  direct  application  of  Newton's 
Law  to  a  free  body  defined  aB  the  material  contained  within  the  chmnber 
and  nozzle  at  time  J,  .  Thms ,  the  force,  Pit),  on  this  free  body  in 
the  tC  direction  is 

cr/  ,  \  -  J—  ‘L('W-U')  <0' 

FU)  -  M; 

/mr clt+At)  -  -w <u, ( -t) 


FU) 


/c  At-»> 


A  Jr 


(4-1) 
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where  Is  the  appropriate  factor  of  proportionality,  m  is  mass  ano  u 
is  fluid  axial  velocit  The  momentum,  mu,  at  time  t  is  approximately 

nvuV)  «  <4“ 

where  f  is  the  fluid  density  and  A  and  I  are  nozzle  dimensions  shown  on 
figure  7,  This  relation  ignores  the  axial  velocity  of  the  fluid 
approaching  the  tube  a.id  the  over  velocity  in  the  "vena  contracts"-  A 
small  constant  correction  factor  could  be  introduced  to  account  for  these 
effects,  if  desired,  with  essentially  no  change  in  the  analytical  results. 

There  are  two  cases  to  consider,  su,  >0  and  'to  <  O  .  The  esse 
for  /vl  >  O  is  treated  first.  Thus,  the  momentum  at  time  t+-$t 

must  account  for  that  portion  of  the  free  body  that  has  left  the  nozzle. 
Since  the  length,  A. I  ,  associated  with  this  portion  is  approximately 

lv(i)  + 


'rft  u  l £ + a.  -  f  I'tALt)  £_]  (momentum  in  tube) 


(momentum  of  portion  outside  nozzle) 


+  terms  of  higher  order; 


'K>C. 


From  (4-1),  therefore, 


Fit) 


p/3  (  A  l 

-2  <  t  rUij 

l  <u 


s\.  \ 


(4-3) 


An  integration  of  forces  acting  on  the  free  body  in  the  +u  direction 
shows  that 


FU)  =  ^  i  ^>0)  <*"*> 

where  p-p^  is  the  difference  in  static  pressure  between  the  motor  chamber 
interior  a.id  the  surrounding  fluid.  F,  is  the  shear  force  on  the  fluid 
passing  through  the  nozzle  which  can  be  expressed  in  terms  of  the  dynamic 
pressure,  o  -  1  f  a.s 

*  ‘  J< 
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£ 


«•*  ;  'U.  >  o 


(4-5) 


whare  CD 

nozzle. 

giving. 


la  the  drag  coefficient  and  S  Is  the  Inner  surface  area  of  the 

Since  a  •  -21  IT  equations,  (4-3)  and  (4-4)  nay  be  equated 
4 


t 


71 


Jco 

g.  db 


■t  £2-  [/  +ZCD  A\  ;  4C  >o 
Jc  1  o  J 


(4-6) 


The  liquid  la  considered  to  be  incormresplble  and,  therefore  any 
Increase  In  air  volume  must  cause  a  corresponding  volume  of  fluid  to  leave 
the  chamber ,  or 


gy  - 


(4-7) 


In  addition,  a  relation  between  the  chamber  pressure,  p,  and  the  air  vol«e, 
V,  may  be  written  as  pVn  ■  a  constant,  where  I  $  m  $  l-  +  .  'yi-t-0 

for  isothermal  and  mm  1.4  for  reversible  adiabatic  (lsentropl**)  processes. 
Prom  this  relation, 

M  v  At  V  <4 

Dlfferantlatlng  (4-6),  while  treating  as  a  constant,  and  substituting 
(4-8)  ylalds 


_  7  to 

$<■  **  *■ 


(4-9) 


Tha  case  for  'ic  <0  is  considered  next.  The  free  body  under 
consideration  is  again  the  tluid  within  the  chamber  and  nozzle  at  time  t. 

The  portion  of  the  nozzle  occupied  by  the  free  body  at  time  t  -r  A  t 
has  a  length  of  approximately 


J  f  (40  +  J  At)  4  /  =  jP  +  stt.  d/ 
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Since  Al.  <  0  ,  this  length  will  be  less  than  /  .  At  time  t+  4  t  the 

momentum  within  the  nozzle  is 


/mutt  +-4  ■-  a i  J 

This  expression  is  used  for  the  momentum  of  the  free  body  under  considera¬ 
tion  and  neglects  the  momentum  of  the  mass  in  transit  between  the  left 
hand  end  of  the  nozzle  end  the  inner  chamber  wall.  Thus,  the  total  force 
in  the  +u  direction  is  again 

F  „  _L_  <L(smno)  -  —  Ltryt 

<£t  3<-  a t—>c  4 £ 


(4-10) 


which  is  identical  with  (4-3) . 

The  evaluation  of  the  force' in  terms  of  pressures  and  shear  force, 
however,  is  somewtut  more  complicated  than  for  the  case  /ic  >  O  . 


Figure  8.  Flow  Into  Chamber 


For  't'C  <  0  the  jet  impact  on  the  inner  wall,  as  shown  in  figure  8, 
causes  a  local  impact  over-pressure  of  an  amount  £p  above  the  chamber 
pressure  p.  Since  there  is  no  net  axial  force  developed  by  the  motor 
when  'u,  <  0  end  -o  ,  i.o.,  steady  state,  the  integral  of 

the  over-pressure  may  be  shown  to  be  equal  to  the  momentum  flow  divided  by 
8c »  or 


ctfi 


cf*' 


-  o 
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The  maes  of  the  fluid  moving  axially  between  the  nozzle  and  the  inner 
chamber  if  small  and  the  jet  transit  time  within  the  chmober  ie  considered 
negligible.  Consequently,  this  expree&ion  is  used  for  the  Integrated  o/er- 
pressure,  even  in  caie  Thus,  letting  represent  the  pressure 

within  tha  nozzle  at  the  edge  of  the  free  body  at  time  t+4t  (see  figure  8) 
an  expression  for  F  in  the  +u  direction  is 


7  ~ 


L?rM 


Hr*}/ 

t  fif'**' 

3+ 


■io 


The  vena  contracts  at  the  inlet  causes  an  inlet  pressure  loss,  and  from  anv 
text  on  hydraulics  it  is  found  that 


t-  -•  '  7T 

or 

/-  —  C  4<?  ~  u  ^ °  -I  ^  ^  ^  ®  ■  (4-11) 

Thus  from  (4-10)  and  (4-11) 


f  -ft*** 


tltO 

<** 


/  v 


zca±] 


'to 


**  <0  (4-12) 


Differentiation  of  (4-12)  and  using  (4-8)  , 


«  ,  * 
a 

Jc  aLt  * 


riL2t-  j  /  +2Cn  1  cLUj  .  mpn-U,  _ 

Jc  L  *  ~v 


Equation*  (4-9)  and  (4-13)  may  be  combined  and  written  as 


(4-13) 


4-  [  t  +  2c.  (Jjl  +  nfflic  c  (4-i4) 

x  Sc  1  &  y 

without  restriction  on  tha  sign  of  u.  Thelulfsctor  in  the  second  ter* 
causes  the  system  to  be  nonlinear  and  not  subject  a  simple  linear  analysis 
of  a  second  order  system.  Note  aleo  that  a  change  in  sign  of  u  rssults 
in  ths  change  in  sign  of  all  terms  on  the  left  hand  side  of  (4-14).  This 
is  to  be  expected  since  the  nozzle  is  eynmetrlc  and  a  reverse  of  th)  pressure 
driving  force  must  reverse  the  flow. 
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( 

An  analogous  relation  to  (4-14)  for  a  series  R,  1,  C  electric  circuit  is 


( 


l  4^-  +  i  --  —■ 

M  *  c  ^ 


(4-15) 


where 

L  ■  inductance,  henrys 
1  "  current,  amperes 

R  ■  resistance,  ohms 
C  ■  capacitance,  farads 
e  ■  electromotive  force,  volts 

The  FPJM  can,  therefore,  be  considered  to  be  the  analog  of  the  nonlinear 
equivalent  circuit  shown  In  Figure  9. 


efr) 


Figure  9.  FPJM  Equivalent  Circuit 
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4.2  ileaonant  Frequency 

While  an  exact  analysis  of  the  nonlinear  system  of  equation  (4-14) 
or  figure  9  is  rather  involved,  the  resonant  frequency,  fc  ,  for  the 
limiting  case  of  R  ■  0  may  be  easily  obtained  as 


~  27T  fZc 


'TfIE*L 

r  irfjrr 


D 


(4-16) 


Using  gc  “  32i/ ^  ■  '2  ~  3g L 


f  -  &Z.  ¥  f  .  — 

J  L  ft'  77 


l _  f/*\  I 

A  <?  3 


O.  e>3c  /  /£*»/■, 


1.4,  for  an  isentropic  process 


p  -  14.7  lbf/in2, 
This  expression  becomes 


/c 


6L-t& 


Z'a' lCt. 


with  units: 


D  (inches)  ,  j[  (inches)  ,  V  (cubic  inches)  and  /(hertz  ■  cycles  per 
second) . 

4 . 3  Approximate  Steady-State  Solution  for  Motor  Discharge  Velocity 

An  approximate  solution  for  the  second  order  nonlinear  differential 
equation,  4-14,  can  be  obtained  by  using  the  describing  function  as  discussed 
by  Kochenburger ,  ref.  2.  With  this  technique  the  input-output  relation 
for  the  nonlinearity  may  be  approximated  as  shown  in  figure  10. 


Tnpu  t 

tfeu  />  f  y 

-u(i)  -  lulu-  '/ 

f  /lo{  ctu> 

fv  -- 

"  /*- 

j/6 

Ou  tptf  f 


itoim 

3?r 


CC5u.'t  - 


4JU\_Au, 

3T  d* 
- 


Figure  10. 


Describing  Function  Approximation  for Nonlinearity 

<fb 
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Using  this  approach,  all  harmonics  of  the  output  which  are  of  higher 
order  than  the  first  are  neglected.  If  the  pressure,  ^  ,  contains  a 
sinusoidal  variation  about  an  average  value,  the  fundamental  component 
of  the  nozzle  discharge  velocity  will  vary  at  the  same  frequency,  uJ/i'iT. 
Thus,  representing  the  amplitude  of  the  velocity  variation  by  J£Tj  , 

r,Aj  -  l  \  St*  nit 


And  equation  (4-14)  may  be  used  to  find  the  required.  A 

particularly  convenient  method  of  solution  uses  phasors  and  in  this  notation, 
along  with  the  describing  function  approximation,  equation  (4-14)  has  the 
steady-state  solution 


7 


7 


or 


U  - 


-  ? 


SflOl  [t*Coi&\\ 

-A-  -A 

l  f**3 

J 

/n  -p  d  \ 


(4-18) 


where  J  ■  /- /  and  upper  case  U  and  P  are  velocity  and  pressure  phasors. 

The  corresponding  phasor  equation  for  the  analogous  electric  circuit  is 


I  - 


(4-19) 


The  describing  function  approach,  therefore,  allows  the  identification  of  a. 
simpler  expression  for 


(4-20) 


Mote  that  the  equivalent  L  and  C  are  unchanged  from  those  identified  In 
figure  9.  From  (4-20^  it  is  seen  that  the  equivalent  resistance  is 
proportional  to  the  amplitude  of  U.  Thus,  the  circuit  selectivity  or  Q, 
where 


Q  = 


LOr  L 

~~fT 


l 

S\V\ [l  +  cj£)\ 


(4-21) 
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Evidently,  for  a  fixed  natural  frequency,  ,  the  Q  of  the  motor  varlc-i 
•irectly  with  £  (for  A.C0L.  <<  £.  )  and  inversely  with  the  velocity 

amplitude.  D 


At  resonance  u>L  ■  /wc  and  the  current  from  (4-19)  depends  on^y  upon 
the  circuit  R.  Analogously,  the  motor  the  Jet  velocity  at  resonance,  UQ, 
depends  only  on  Req  and  from  (4-20) 


If, 


-  3  g<f JPj 
Sslvil'+JCuA.] 


or 


(4-22) 


Thus,  the  describing  function  approach  yields  simple  expressions  for  the 
selectivity,  Q,  and  Jet  velocity  of  a  FPJM  at  resonance. 


4.4  Motor  Thrust 


The  average  net  thrust  of  the  motor  at  any  general  frequency,  f,  is 
equal  to  the  average  value  of  the  momentum  flew  leaving  the  device  for  a 
period,  T  ■  1/f.  Thus 


t*  T 


Pavg  -  rp 


AO(0  >  O. 


(4-23) 


or  considering  only  the  fundamental  component  of  u(t), 

I  VI* 

favg  -  *  (4-24) 

.P  + 

At  resonance,  using  (4-22) 


Psvg.c 


S’TTA  ip) 


(4-25) 


The  thrust  at  frequencies  other  than  resonance  can  be  obtained  by  solving 
(4-18)  for  trj|\  First  N  is  defined  as 


(4-26) 
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Thei,  from  (4-18) 

I  Pi" 

,U|  '  VNc/l"  r  [»L-  f 

or 

A/i  [  icrri* +  [  U)L~  ill  /UI*  -  ipi‘--o 


(4-27) 


(4-28) 


Equation  (4-28)  may  be  solved  for  Itfl*"  from  the  quadratic  formula  using  the 
+  sign  on  the  radical  term  since  Ji/i’" must  always  be  positive.  The  average 
thrust  is  then  given  by  (4-24).  A  program  (Appendix  A)  was  written  for  the 
Wrng  700  which  provides  rapid  solutions  for  thrust  as  a  fuuction  of  motor 
volume ,  V,  nozzle  diameter,  D,  nozzle  length,  X  ,  and  amplitude  of  the 
sinusoidal  external  pressure  variation,  F.  The  program  is  based  on  water 
with  a  density  of  62.4  lbm/ft^,  C&  ■  0.01,  an  average  motor  chamber 
pressure  of  14.7  lbf/in2and  n  ■  1.4  for  isentropic  changes  in  air  volume. 
Figures  11  through  14  show  the  variation  of  thrust  as  a  function  of 
frequency  with  different  combinations  of  parameters  held  constant  such 
that  the  undamped  natural  frequency,  fQ  •  10  cycles/sec.  Figure  15 
gives  the  calculated  performance  of  the  port  and  starboard  motors  of  the 
submarine  model  shown  in  figure  6. 
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Figure  11.  Jet  Ho  tor  Thrust  with  Variations  in  Nozzle 
Tube  Length  and  Chamber  Volume 


CALCULATED  THRUST 


Figure  12.  Jet  Motor  THrust  with  Variations  in  Pressure  Amplitude 


[otpr  Thrust  with  Variations  in  Nozzle  Tube 
Diameter  and  Chamber  Volume 


Motor  Thrust  with  Variations  in  Nozzle  Tube 
Length  and  Diameter 


Chamber 
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SECTION  V 


DEPTH  STABILITY  OF  CARTESIAN  DIVER 


A  simple  Cartesian  Diver  has  an  equilibrium  point  at  a  depth  determined 
by  the  amount  of  trapped  air  and  the  mass  of  the  device.  This  equilibrium 
point,  however,  is  unstable  whenever  it  exists  between  the  upper  and  lower 
surfaces  of  the  tank  containing  the  diver.  Thus,  if  the  diver  position  is 
perturbed  either  upward  cr  downward  from  such  a  point,  a  force  is  developed 
in  a  direction  to  increase  the  perturbation.  If,  for  example,  the  device 
is  forced  below  the  equilibrium  point,  the  hydrostatic  pressure  on  the 
trapped  air  becomes  higher,  the  air  volume  becomes  smaller,  buoyancy 
decreases  and  the  device  tends  to  sink  further,  provided  no  correction  force 
is  applied.  Such  a  correction  force  can  be  furnished  by  reducing  the 
pressure  of  the  surrounding  fluid  as  in  the  case  of  the  familiar  toy. 
Alternatively,  a  corrective  vertical  jet  force  can  be  controlled  by  a 
suitable  sensing  element  within  the  diver.  Ona  such  possible  control 
device  is  shown  in  figure  16  and  described  in  the  following  section. 

5 . 1  Description  of  a  Depth  Control  System  (DCS) 


Figure  16.  Depth  Control  System  (DCS)  for  Cartesian  Diver 
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Operation  of  the  depth  control  system  (DCS)  of  figure  16  is  accomplished 
by  throttling  the  nozzle  of  a  fluidic  pulse  jet  motor  (FPJM)  with  a  flexible 
diaphragm.  Assume  that  the  fluid  surrounding  the  diver  is  pulssting  sbout 
some  fixed  average  value  and  that  the  diver  is  at  the  equilibrium  point. 

This  insures  that  theve  is  an  equality  between  the  downward  gravity  force 
on  the  diver  and  the  sum  of  the  upward  forces  due  to  both  the  buoyance  of 
the  air  space  and  the  FPJM  average  thrust.  When  the  diver  position  is 
perturbed  downward  there  is  again  an  increase  in  average  pressure  on  the 
air  space  and  a  decrease  in  buoyancy  due  to  the  resulting  reduction  in 
Bir  volume.  With  the  DCS  in  operation,  however,  the  decrease  in  air 
volume  results  in  an  opening  of  the  pssssge  between  the  top  of  the  FPJM 
nozzle  tube  and  the  flexible  diaphragm  which  increasis  the  vertical  jet 
thrust.  When  this  incresse  in  Jet  thrust  is  greater  thsn  the  loss  of 
buoyance,  there  will  be  a  net  corrective  force  tending  to  eliminate  the 
initial  downward  position  perturbstion.  It  is  also  easy  to  see  thst  a 
downward  corrective  force  will  arise  if  the  initial  position  perturbation 
from  equilibrium  is  un«,crd. 

Quantitative  design  requirements  to  insure  thst  a  DCS  will  provide  a 
stable  equilibrium  are  developed  in  the  next  section  of  the  report. 

5.2  Analysis  of  Depth  Control  System 

The  DCS  for  the  Cartesian  Diver  shown  in  figure  16  may  be  analyzed  by 
first  considering  the  region  between  the  nozzle  tube  and  the  flexible 
diaphragm.  Figure  17  shows  the  diaphragm  and  the  upper  end  of  the  nozzle 
separated  by  a  distance  x.  It  is  assumed  that  the  pressure  in  the  tank 
containing  the  diver  is  pulsating  sinusoidally  with  a  fixed  frequency  and 
amplitude  about  sn  average  value  and  that  equilibrium  is  obtained  when 
the  time-average  value  of  x,  designated  by  x,  is  equal  to  Xq. 


DfapMaxp'rrt 

trv 


X 


Figure  17. 

Detail  of  Region  Betweeu  Nozzle  Tube 
snd  Diaphragm 
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It  Is  also  assumed  that  the  upward  jet  thrust  on  the  diver  caused  by  the 
pulsating  flow  leaving  the  nozzle  monotonically  increases  with  increasing  x. 
When  x  “  o,  no  flow  can  pass  through  the  nozzle  and  thus  no  upward  thrust 
results.  As  x  Increases  from  zero,  the  upward  thrust  increases  until  until 
x  reaches  a  value  comparable  to  the  nozzle  inner  diameter  when  the  diaphragm 
provides  essentially  no  nozzle  throttling  and  the  FPJM  thrust  approximates 
the  value,  calculated  by  the  method  described  in  Section  4. A.  Although 

an  exact  formulation  of  the  relation  between  the  opening,  X,  and  the  upward 
motor  thrust,  Fm  is  not  available,  it  is  expected  to  appear  as  shown  in 
figure  18, 


Figure  18.  FPJM  Thrust,  Fm,  vs  Nozzle  Tube  Opening  x. 
C Dm  is  the  nozzle  inner  diameter,  Xq  is  the  opening  at 
equilibrium  and  Fm0  is  the  corresponding  thrust). 


At  the  equilibrium  point  the  total  net  time-average  vertical  force  in 
the  +y  direction  (upward)  must  be  zero.  Three  such  forces  act  on  the  diver, 
the  downward  force  of  gravity,  the  upward  fluidic  pulse  jet  motor  (FPJM) 
thrust  and  the  upward  buoyant  force  equal  to  the  weight  of  the  displaced 
liquid.  The  latter  two  forces  have  both,  (1)  a  periodic  component  with  the 
same  frequency  as  the  external  pressure  variation  and  a  mean  value  of  zero 
and.  (2)  a  longer  term  "average  value"  about  which  the  periodic  component 
varies.  The  periodic  component  of  these  forces  causes  a  small  amplitude 
periodic  motion  of  the  diver  to  be  superimposed  on  the  motion  caused  longer 
term  changes  -in  average  force  values.  The  small  periodic  effects  are 
neglected  in  the  following  analysis  and  attention  is  concentrated  on  the 
longer  term  effects  caused  by  changes  in  average  values  of  forces  and  diver 
position,  y. 

Expressions  relating  the  average  values  of  the  vertical  forces  and 
position  of  a  diver  may  be  derived  by  considering  upward  values  as  positive 
Thus  the  force  of  gravity,  Fg,  in 

(3  (>  **/•/> 1 )  . 

F.  (M.f)  *  -  ///*./ 

3  y<  f  j't*  f-y-t  /  if 
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The  buoyant  force,  Fb ,  depends  on  the  total  volume  of  displaced  fluid, 

Vs  +  Va,  where  Vs  is  the  volume  of  the  solid  parts  of  the  diver  and  Va 
is  the  air  volume.  Thus 

jc  v,>  - 

Letting  FOT  represent  the  FPJM  thrust  and  using  the  subscript  "o"  to  indicate 
equilibrium  conditions, 


F  +  fjt  A  ^  ° 

or,  since  g/g c  has  a  numerical  value  of  1  st  sea-level 


(5-1) 


**7  *■  f  ^  ^  r  %  * 


+  F. 


--o . 


(5-2) 


When  the  diver  •'s  perturbed  from  the  equilibria  depth,  yo,  by  an 
amount  $y,  the  gravity  '->rce  is  unchanged.  However,  the  average  FPJM  thrust 
will  change  to 


1 


(5-3) 


A  reasonable  value  for  ^Fsc  may  be  obtained  by  considering  that  the  thrust 

changes  from  zero  to  full  value,  F^,  as  x  changea  from  zero  to  Dn  or  as  an 
approximation , 


'/ftr.'+vf 

o 


(5-4) 


Letting  Dj  be  the  effective  diameter  of  the  diaphragm  shown  in 
figure  5-1,  the  t4me  average  change  in  volume,  dVa,  of  the  airapace  is 


f  (*J) 


(5-  j) 
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From  the  usual  expression  for  the  air  volume 


f 


K 


<?  Cot* 


5 /-a 


c/l't 


where  the  bar  over  ti  .  symbol  p  again  denotes  a  time  average. 
F-om  the  hydrostatic  relation, 

Jf  --  - • 


Thus 


,  _  jfs.'k 

Jv  -'S 


and  from  equation  (5-3) 


-  t 


<4. 


Note  that  equation  (5-8)  defines  the  constant  ■*. , 
and 

> 


-  /r  -  oi-  Xt 


*rro 


Letting  represent  the  buoyant  force 

'  (*&)** 


where 


J 


or 


,c  I  ,  JjA 

^  -  rU 

jr<  ,  -  ?4~-4 -  ■ 


r 


*3 


Ar  JT 


Use  of  relations  (5-6)  and  (5-7)  shows  that 


tft 


x  /  a  - 

j>  ^  3 

/yt 


(5-6) 


(5-7) 


(5-8) 


(5-9) 


(5-10) 
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which  defines  the  quantity  0  ,  Thus 


(5-11) 


Assuming  a  viscous  drag  force  of  an  amount  - 


the  differential 


equation  of  motion  for  small  changes  from  equilibrium  becomes 


/ 


t  *■ 


c\- 


jr 


r 

A- 


J. 


1 


P 

.'•Pc 


o L 


(5-12) 


The  value  for  "m"  used  in  this  relation  must  include  any  liquid  within  the 
diver  that  will  be  accelerated  along  with  the  actual  structure  of  the 
device,  Because  of  (5-1)  the  equilibrium  terms  in  (5-12)  may  be  eliminated 
and  this  relation  rewritten  as 


f  cm  •[_ 

t  S3  /•=•)  r  -u 

■c 

L  ’ 

M  i  J 

(5-13) 


For  stability,  all  coefficients  in  this  expression  must  be  positive.  Thus, 
it  is  necessary  that 

o(  >6 

Or,  from  (5-8)  and  (5-10) 


f 

sm  w 


/> 


'J 


(5-14) 


Relation  (5-14)  states  the  requirement  for  depth  stability  of  a  oiver  with 
a  FPJM  DCS.  As  noted  in  the  analysis  it  is  based  on  an  approximation  for 
the  variation  of  FPJri  Uuust  w  iiu  v  i aphr no  zzle  opening;  It  Rhould  also 
be  observed  that  ,  the  effective  diameter  of.  the  diaphragm,  is  actually 
defined  by  equation  (5-5)  through  the  relation  between  air  space  volume 
Va,  and  diaphragm-nozzle  opening  x.  may,  therefore,  be  difficult  to 
commute  for  a  given  design.  In  view  of  these  facts,  relation  (5-14)  must  be 
considered  more  of  a  guide  indicating  the  important  factors  for  stable  operation 
rather  than  a  go-no-go  design  criterion.  The  relation v  for  example, 
indicates  that  as  the  amplitude  of  pressure  pulsation  increases  from  zero 
and  thrust,  Fog,  ,  increases  a  value  will  be  reached  for  which  depth  stability 
is  first  attained.  Stability  will  then  be  maintained  for  higher  pressure 
amplitudes  at  the  same  pulsation  frequency.  The  relation  also  indicates  the 
desirability  of  a  small  effective  diaphragm  diameter  which  is  equivalent 
to  a  large  diaphragm  movement  for  a  given  change  in  fluid  pressure. 
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SECTION  VI 

SYSTEM  MEASUREMENTS 

This  section  describes  measurements  made  of  tank  pressure  variation, 
motor  thrust  and  submarine  velocity. 

6 . 1  Pressure  Variation 


Tank  pressures  were  measured  using  a  pressure  sensitive  transistor 
(P1TRAN)  made  by  S tow -Labor ator i es ,  Inc.  The  PITRAN  was  mounted  in  the 
tank  hstch  and  calibrated  against  a  water  filled  manometer  as  in  figure  19. 


Vent- 


71 

in  .out  p 

o/Kh  ip 
auoospne^  e 


^  Control 


Vidte'r 

Mmme/e't' 


Figure  19.  P1TRAN  Calibration  Set-up 
The  electrical  circuit  is  shown  in  figure  6-2. 


31 


NAVTRAEOUIPCEN  IH-187 


Sufficient  air  was  retained  in  the  top  of  the  ta\k  so  that  the  pressure 
control  piston  could  be  fully  ir.serted  without  forcing  water  out  the  tank 
vent  and  into  the  manometer.  When  the  piston  was  withdrawn  the  water  in 
the  tank  decreased,  perhaps  2  inches,  as  air  from  the  manometer  tube  entered 
the  tank.  This  caused  a  change  in  hydrostatic  pressure  on  the  PITRAN  of 
less  than  0.1  psi.  The  calibration  curves  shown  on  figure  21  t  however, 
are  based  on  the  pressure  measured  by  the  manometer.  Thus,  it  was  assumed 
that  changes  in  pressure  at  the  tank  vent  were  an  adequate  measure  of 
changes  in  pressure  at  the  PITRAN,  and  the  0.1  psi  change  due  to  tank 
level  variation  was  neglected. 

From  figure  21  ,  the  PITRAN  was  found  to  have  a  calibration  of  3.2 
psi  per  volt  (0.31  volts/psi)  when  ■  1  volt  at  zero  differential 

pressure.  The  factory  calibration  was  1.73  psi  per  volt  (0.58  volts/psi) 
with  Vq£  =  2  volts.  Stresses  introduced  into  the  transducer  by  the  tank 
mounting  hardware  may  have  been  responsible  for  the  observed  loss  in 
sensitivity.  Sufficient  transducer  sensitivity,  however,  was  available 
to  measure  tank  pressure  as  shown  on  figures  22  through  28  which  were 
all  recorded  using  a  CRO  set  at  0.05  volts/cm  on  the  vertical  axis  and 
20  ms/cm  on  the  horizontal  axis.  The  vertical  axis  pressure  calibration 
is  therefore,  3.2  x  0.05  “  0.16  psi/division. 

Maximum  amplitude  of  the  piston  periodic  i.otion  is  possible  only 
when  the  average  position  is  near  mid-stroke.  Near  either  end  of  the 
stroke  periodic  motion  is  limited  by  the  control  system  to  avoid  impact 
of  the  hydraulic  piston  with  the  ends  of  the  cylinder,  see  p.  52  of 
refeience  1.  This  "limiting"  causes  the  difference  between  the  traces 
on  figure  25  and  27  .  The  amplitude  control  setting  was  unchanged  for 
the  two  conditions  and  the  only  difference  was  that  the  average  piston 
control  setting  changed  from  1/2  to  3/4  withdrawal,  i.e.,  from  2.6  inches 
to  3.75  inches  outward  from  full  insertion. 

One  observed  difficulty  with  the  system  was  that  two  different 
amplitudes  of  motion  occurred  with  the  same  amplitude  control  setting. 

For  example,  all  control  settings  for  figures  23  and  24  were  the  same 
but  the  system  exhibited  two  modes  of  operation.  The  change  from  one 
piston  motion  amplitude  to  the  other  was  abrupt  and  one  mode  appeared 
to  be  stable  until  the  change  to  the  other  occurred.  The  same  phenomena 
was  observed  in  the  change  from  figure  25  to  26  .  Neither  mode  of 
operation  teemed  to  occur  more  often  than  the  other. 
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Figure  22  .  Hydraulic  Pump  On.  Average  Piston  Position: 
1/2  Stroke.  Piston  Amplitude:  0. 
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Figure  23.  Average  Piston  Position:  1/2  Stroke.  Piston 
Motion:  to.25  inch.  Pressure  Variation:  0.12 
psj .  f  *  10.6  cycles/sec. 


Figure  24.  Average  Piston  Position:  1/2  Stroke.  Piston 
Motion:  —  0.5  inch.  Pressure  Variation: 

0.21  psi.  f  *  10.6  cycles/sec. 


Figure  23.  Average  Piston  Position:  1/2  Stroke.  Piston 
Motion:  £  0.28  inch.  Pressure  Variation: 

£  0.21  psi.  f  a  8  cycles/sec. 
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Figure  26.  Average  Piston  Position:  1/2  Stroke.  Piston 

Motion:  —0.5t  inch.  Pressure  Variation:  £  0.37 
psi.  f  *  8  cyclee/sec. 


Figure  27.  Average  Piston  position:  3/4  Stroke  (out). 

piston  Motion:  —0.16  inch.  Pressure  Variation 
—  0.09  psi.  f  *  8  cycles/sec. 
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Figure  26  ,  Average  Piston  Position:  3/4  Stroke  (out). 

Piston  Motion:  i  0.34  inch.  Pressure  Variation 
—  0.22  psi.  f  =  8  cycles/sec. 
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6.2  Motor  ThruBt  Measurement 

Figure  29  18  a  photograph  of  a  model  pulse  jet  motor  which  was  used 

to  measure  thrust.  The  motor  had  an  air  space  of  0.3  In3  with  a  nozzle 
tube  0,425  In.  long  and  1/16  In.  Internal  diameter.  The  model  was  attached 
to  a  stainless  steel  ball  chain  such  as  Is  often  uaed  for  key  chains.  The 
chain  weighed  4.8  grams  In  air  and  3.7  grams  in  water.  It  was  10  In. 
long  and  was  made  up  or  62  balls  and  links.  Thus,  in  water,  it  required 
0,13xl0“3  lbf  to  raise  one  Hell 


Figure  29,  Thrust  Model 

from  the  tank  floor.  Measurements  showed  that  at  11  hAttz  a  pressure 
pulsation  of  0.32  psl  amplitude  raised  2  balls,  indicating  that  the  motor 
was  developing  0.26xl0-3  lbf.  The  thrust  program  Included  In  the  appendix 
yields  a  resonant  frequency  of  11.6  hertz  and  a  thrust  of  0.25xl0"3  lbf 
at  the  conditions  used  for  the  measurement.  The  agreement  between 
experiment  and  theory  was  surprisingly  good. 

6.3  Model  Velocity 


A  computation  of  expected  model  velocity  was  made  based  on  the  motor 
dimensions  given  on  figure  6.  The  model  was  estimated  to  have  a  surface 
area  of  24  in^  and  a  42°  angle  between  jet  and  model  axes.  The  reault 
was  a  prediction  of  1.8  ft/sec,  or  22  in/sec.  At  this  velocity,  the 
Reynolds  Number  was  9.7xlo\  the  drag  coefficient  from  the  Blaslus 
formula,  see  Eshbach  -  reference  3,  page  6-37,  was  0.00443,  the  dynamic 
pressure  was  3.14  lbf/ ft2,  and  these  conditions  yield  a  computed  drag  of 
2. 3x10 *3  lbf. 


At  a  frequency  of  11.4  hertz  and  a  pressure  amplitude  of  0.37  lbf/in 
ths  thrusts  from  the  port  and  starboard  motors  are  equal  at  a  value  of 
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1.56x10"^  lbf.  Thus,  since  cos  42°  ■  0.74,  the  total  available  thrust  is 
computed  to  be  2.3x10"^  lbf. 

Physical  measurements  of  model  velocity  proved  to  be  impossible 
because  of  the  limited  tank  dimensions  and  the  difficulty  of  maneuvering 
the  model.  The  highest  velocity  actually  observed  by  timing  the  model 
between  two  lines  spaced  one  foot  apart  on  the  tank  was  1/2  ft/sec.  The 
model  was  accelerating  and  might  have  attained  a  steady-state  velocity 
in  the  neighborhood  of  the  computed  22  in/sec. 
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SECTION  VII 

CONCLUSIONS  &  SUGGESTIONS  FOR  FUTURE  INVESTIGATIONS 

A  theory  of  operation  of  a  Fluidic  Pulse  Jet  Motor  (FPJM)  was  developed 
which  allows  computation  of  trotor  thrust  by  means  of  a  computer  program. 
Excellent  agreement  was  established  between  computed  and  measured  thrust 
values . 

A  proposed  design  of  a  depth  control  syatem  (DCS)  was  devised  which 
should  be  applicable  to  any  model  using  the  Cartesian  diver  principle. 

Depth  stability  relations  were  developed  and  used  to  obtain  a  design 
criterion  for  a  DCS. 

Based  on  experience  with  the  thrust  measuring  model,  figure  29,  ,  a  new 
model,  figure  30  ,  was  designed  which  can  be  used  to  measure  thrust  with 
two  alternate  nozzles.  The  new  model  can  also  be  used  to  verify  the  depth 
control  system  principle  by  adding  a  flexible  rubber  diaphragm  above  the 
upper  end  of  the  nozzle  as  shown  in  the  figure. 

It  is  recommended  that  the  new  design  be  built  snd  used  for  thruat 
and  depth  stability  measurements.  Another  worthwhile  project  will  be  to 
design  and  build  a  horizontally  maneuverable  model  using  a  suitable  depth 
control  system.  Following  a  successful  demonstration  of  the  latter  device, 
a  practical  training  device  csn  be  undertaken. 


Figure 
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APPENDIX  A 


PROGRAM  FOR  FLUIDIC  POi.SE  JET  MOTOR  THRUST 


PROGRAM  TITU. 

NUMBER 

N.A. 

FLUIDIC  PULSF.  JET 

MOTOR  (FPJM)  THRUST 

PROGRAMMED  BY 

PROGRAM  ABSTRACT 

H.  C.  TOWLE 

COMPUTES  THRUST, 
SELECTIVITY  (0) 

NATURAL  FREQUENCY  6 

DATE  j 

JULY  1 9  /  ;> 

1 

1 BLOCKS 

NO.  OF  STEPS  j 

NO. 

NO.  OF  STEPS 

DATA  REGISTERS 

MARK  USED 

225 

1 

22: 

00-15 

Mark  0,i 

VERIFY  NUMBER 

2391  s 

SET  P.C.  J  PAT 

000  l  1  of  7  J 

i 

i 

Load  ctifc  progr.-' 

tape.  Depress  LEV 

TND;  Depruvo  RUN 

TAPE  ’  .  '  , 

PRIME;  LOAD  PROGRAM;  PRIME;  LOAD  PROGRAM  (#2  program 

on  tape)  .  ■ 

i 

2. 

PRIME;  VERIFY 

PROGRAM;  x  register  should  read  2391. 

3. 

SEARCH;  0.  Computes,  displays  and  stores 

j 

/7 gc  (Ibm/ in3) . 

(  y  :  -f-  .935424078100  -04 

f  -fluid  density,  lbm/in3 

(x  :  \  .935424078100  -04 

Sc  *32.17 

•f  12,  dimensionless 

1 

\ 

A  =  TFD2 

4 

— 

1 

A  -  nozzle  tube  area  in^. 

D  -  nozzle  tube  diameter,  in 

■ 

. 

%  -  nozzle  tube  length,  in. 

i 

4. 

Key  l;  ♦  ;  D; 

GO. 

EXAMPLE 

0.425;  ♦  ; 

i 

.0625;  GO 

Read  y  :  A 
x  :  4 

y:  +  .3067961575/5  -02 
x:  4  4.00. . , 
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PROGRAM  TITLE 

FLUIDIC  PULSE  JET  MOTOR  (FFJM)  THRUST _ PAGE  2  of  9 

l  «  £i_ 

8c 

C  -  V 
nAp 

N  *  iiL  f  1  4-  2Cp  1,1  =  R 

3^c  l  D  J  |U| 

V  _ 1 

2  tt  rrr 

C  -  equivalent  capacitance,  in^/lbi, 

Op  -  drag  coefficient, assumed  =  0.01,  dimensionless. 

f0  -  undamped  natural  frequency,  cycles/aec. 

L  -  equivalent  inductance,  lbm/in^* 

N  -  as  defined  above,  lbm/in^, 
p  -  average  motor  chamber  preasure,  lbf/in^. 

R  -  equivalent  reaistance,  lbm/in^  -  aec. 

I U i  -  amplitude  of  sinusoidal  jet  velocity,  ft/sec. 


EXAMPLE  (continued) 

Key  V;  GO 

Key  0.393:  GO 

n  .  .  j  ..  .  c 

ivcou  y  .  Lq 

Read  y:  4-  10.1174895006 

x:  1/ f0 

x:  .988387484799  -01 

005:  L 

005:  4-  .397555233192  -04 

006:  C 

006:  4-  6.22439670150 

007:  N 

007:  -f  .901993334297  -04 

F0  •  3*T  AP _  , 

32^  +  2CD^j 

F0  -  motor  average  thrust  at  resonance,  lbf. 

P  -  amplitude  of  sinusoidal  change  of  preasure 
in  surrounding  fluid,  lbf/in^. 
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F  =  lk_  |'J|2 

Q  »  0>oL  • 

R 

a  s  N2 

b  •  j«L  ■■  1/fcjcj  2 

c  =  -  P2 

F  -  average  thrust  of  motor,  lbf 

f  -  frequency  of  sinusoidal  change 

of  pressure  in  surrounding  fl«3d,  cycles/sec. 
Q  -  selectivity,  dimensionless 


'0  =  2  7Tf  ;  eOQ  =  2  TTf0 


1 

£, AAIU LZ  ^LUUL  ItlUCu) 

7.  Key  f;  GO 

Key  5;  GO 

Read  y:  F 

Read  y:  +.654457482480  -03 

x:  f 

x:  +5.000. . . 

013:  1  Ut 2 

013:  +9121,85087551 

014:  q 

014;  +.293361343431 

EXAMPLE  (continued) 

8,9,...  Repeat 

Key  20  ;  GO 

Step  7  for 

all  desired 

Read  y  :  +.658957612802  -03 

values  of  f. 

x  :  +20.00... 

013  :  +9184.57384655 

L. - 

0i4  :  +.292357920840 
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Slip 

Key 

(lode 

Comment 

U'litiSmnnsnsHHH  1 

ConmnU 

160 

t 

,0404 

-c  ' 

-c 

V  1 

HE  DIR 

Hi  2 

007 

000/ 

N 

■ 

1..3 

mmm 

Q/1,3 

a  •  N~ 

it 

m 

01,02 

-ac 

DB 

4 

07(1'. 

7. 

|  1  * .  Ci 

X 

0..07 

-  %*J  c 

r.-7 

Old  1 

— '«i*C 

H3 

Dll! 

O'lOO 

m 

.  . 

■mi 

-4ac 

1 1  0 

KE  DIK 

O',  06 

17  1 

012 

0012 

b^-^ac 

1/  2 

0.  12 

WbJ-4ot 

MPi 

Ufil 

1 

012 

012 

't 

RE  DIR 

0,05 

1/  6 

st  g 

VB" 

17  7 

i n 

0/13 

.  b 

■1 

178 

-  DIR 

0401 

17  9 

012 

0012 

-h  +^, 

iu  o 

„ 

4. 

0/  02 

** 

1 

it)  i 

t 

0,.07< 

2 

13  2 

RE  DIR 

0  <06 

13  3 

007 

CC07 

N 

13  4 

X2 

1WU 

A 

13  5 

X 

2a 

la  6 

_ 5' 

0..06 

hsm 

13  7 

V 

.0415 

18  8 

012 

0012 

wr*r 

13  9 

• 

• 

0:03 

Ivl2 

g- *-*-•»•* 

■mm 

19  0 

STORE  Y 

0414 

un 

013 

mi 

■am 

19  2 

RECAIL  Y 

0415 

COMPUTE 

19  3 

003 

0003 

Q 

*n  ... 

19  4 

2 

0702 

2 

19  B 

X 

2f„ 

196 

?r 

H^-il 

7T 

m 

X 

B . 

\'p  — 

fc>o 

RE  IiIR 

0404 

r 

LLLL 

00  > 

00, 1. 

t. 

50 


NAVTRAEOUIPCEN  IH-187 


700  PROGRAM  CODES 


CODE 

KEY 

CODE 

KEY 

0400 

+ DIRECT 

0600 

+ 

0401 

-  DIRECT 

0601 

- 

0402 

x DIRECT 

0602 

X 

0403 

-r  DIRECT 

0603 

0404 

STORE  DIRECT 

0604 

t 

0405 

RECALL  DIRECT 

0605 

0406 

OdRECT 

0606 

0 

0407 

SEARCH 

0607 

IXI 

0408 

MARK 

0608 

INTEGER  X 

0409 

GROUP  1 

0609 

it 

0410 

GROUP  2 

0610 

Lo<3.oX 

0411 

WRITE 

0611 

Log  X 

0412 

WRITE  ALPHA 

0612 

Vx 

0413 

END  ALPHA 

0613 

10x 

0414 

STOREY* 

0614 

X 

C 

0415 

RECALL  Y* 

0615 

2. 

X 

0500 

+  INOIR 

0700 

0 

0501 

-  INDIR 

0701 

1 

0502 

xINDIR 

0702 

2 

0503 

■f  INDIR 

0703 

3 

0504 

STORE  INDIR 

0704 

4 

0505 

RECALL  INDIR 

0705 

5 

0506 

OlNDIR 

0706 

6 

0507 

SKIP  if  Y>X 

0707 

7 

0508 

SKIP  if  Y  <  X 

0708 

8 

0509 

SKIP  if  Y  3  X 

0709 

9 

0510 

SKIP  if  ERROR 

0710 

SET  EXP 

0511 

RETURN 

0711 

CHANGE  SIGN 

i  0512 

ENO  PROG 

0712 

DECIMAL  POINT 

0513 

LOAO PROG 

0713 

XJ 

0514 

GO 

0714 

RECALL  RESIDUE 

0515 

STOP 

0715 

CLEAR  X 

•ENTERED  BY  TOGGLE  SWITCHES 


